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Abstract 

Using the trace anomaly and low energy relations, as well as the Witten- 
Veneziano formula for the mass of the rj meson, we have developed a formal- 
ism which makes it possible to express the gluon condensate, the topological 
susceptibility and the mass of the rf meson as a functions of the truly nonper- 
turbative vacuum energy density, which is one of the most important charac- 
terists of the QCD true ground state. It was directly applied to the numerical 
evaluation of the chiral QCD vacuum topological structure within its quan- 
tum zero momentum modes enhancement model. A rather good agreement 
with the phenomenological and experimental values of the above-mentioned 
quantities has been achieved. With the help of the Witten-Veneziano for- 
mula, we derived an absolute lower bounds for the pion decay constant and 
the mass of the rf meson in the chiral limit. By introducing the most general 
parametrization of the gluon condensate, we also proposed how the correct 
Nf (number of flavors) dependence of its phenomenologically extracted value 
could be restored. 

PACS numbers: 11.15.Tk, 12.38.Lg 
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I. INTRODUCTION 



It is well known that one of the most important aspects of the famous U (1) problem [1,2] 
is the large mass of the 7/ meson. It does not vanish in the chiral limit, so the rj meson is 
not the Nambu-Goldstone (NG) boson. In Ref. [3] (see also Ref. [4]) by using the large N c 
limit technique the expression for the mass of the 7/ meson was derived, namely 

m% = — ^Xt + A, (1.1) 

where A = 2m 2 K — m 2 , Nf is the number of light quarks and F n is the pion decay constant. 
However, the important quantity which enters this formula is the topological density operator 
(topological susceptibility), Xt (for definition see section 3). In the chiral limit it is screened 
that is why it is defined for Yang-Mills (YM) fields, i.e., for pure gluodynamics (Nf = 0). 
It is one of the main characteristics of the QCD nonperturbative vacuum where it measures 
the fluctuations of the topological charge. 

The precise validity of the Witten-Veneziano (WV) formula (1.1) is, of course, not com- 
pletely clear because of its origin. Nevertheless, let us regard it as exact for simplicity (in any 
case we have nothing better than Eq. (1.1). However, there are phenomenological reasons 
[5,6] as well as some lattice indications [7] to believe that QCD is close to SU (00)). Using 
now experimental values of all physical quantities entering this formula, one obtains that 
the phenomenological ("experimental") value of the topological susceptibility is 

^phen = 00X058 GeV 4 = (180.36 MeV) A = 0.1377 GeV/fm 3 . (1.2) 

In the chiral limit A = since K ± and 77 particles are NG excitations. It is worth noting 
further that neither the mass of the 77' meson nor the pion decay constant in the chiral limit 
cannot exeed their experimental values. So the WV formula (1.1) provides an absolute lower 
bounds for the pion decay constant and the mass of the 77' meson in the chiral limit, namely 

854 < m° v , < 957.77 (MeV), (1.3) 

and 

83.2 < F° < 93.3 (MeV), (1.4) 

respectively. They should be compared with their experimental values (upper bounds in 
the previous expressions). Let us note that the chiral perturbation theory value of the pion 
decay constant in the chiral limit, F° = (88.3 ± 1.1) MeV [8], obviously satisfies these 
bounds, Eq. (1.4). Recent lattice result [9] (see also brief review [7]) for the mass of the 77' 
meson in the continuum chiral limit is mj, = 863(86) MeV , satisfying thus bounds (1.3). 

One can conclude in that the mass of the 77' meson remains large even in the chiral limit, 
which is real problem indeed. Thus the large mass of the 77' meson in the chiral limit is due to 
the phenomenological value of the topological susceptibility. In other words, it is clear that 
through the topological susceptibility (i.e., via the WV formula (1.1)) the large mass of the 7/ 
meson even in the chiral limit is determined by the topological properties of the QCD ground 
state, its nonperturbative vacuum. It has a very rich dynamical and topological structure 
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[10-12]. It is a very complicated medium and its dynamical and topological complexity means 
that its structure can be organized at various levels (quantum, classical). It can contain many 
different components and ingredients which may contribute to the truly nonperturbative 
vacuum energy density (VED). It is well known that the VED in general is badly divergent 
[13], however the truly nonperturbative VED is finite, automatically negative and it has no 
imaginary part (stable vacuum). For gauge- invariant definition and concrete examples see 
recent papers [14,15] (for brief description see also sections 2 and 3). Precisely this quantity 
is one of the main characteristics of the QCD true ground state and precisely it is related to 
the nonperturbative gluon condensate via the trace anomaly relation [16] (section 3) as well 
as to the above-mentioned topological susceptibility via the low energy "theorem" (relation) 
derived by Novikov, Schifman, Vanshtein and Zakharov (NSVZ) a long time ego [17] and 
rederived quite recently by Halperin and Zhitnitsky (HZ) [18] (section 4). Let us remind that 
the truly nonperturbative VED is nothing else but the bag constant apart from the sign, by 
definition [13,14,19]. It is much more general quantity than the string tension because it is 
relevant for light quarks as well. 

II. ZMME QUANTUM MODEL 

Many models of the QCD vacuum involve some extra classical color field configurations 
(such as randomly oriented domains of constant color magnetic fields, background gauge 
fields, averaged over spin and color, stochastic colored background fields, etc.) and ingre- 
dients such as color-magnetic and Abelian-projected monopoles (see Refs. [10,11,20] and 
references therein). The relevance of center vortices for QCD vacuum by both lattice [21] 
and analytical methods [22] was recently investigated as well. However, the most elabo- 
rated classical models are the random and interacting instanton liquid models (RILM and 
IILM) of the QCD vacuum. They are based on the existence of the topologically nontrivial, 
instanton-type fluctuations of gluon fields there, which are nonperturbative, weak coupling 
limit solutions to the classical equations of motion in Euclidean space [23] (and references 
therein). That is instantons may be qualitatively responsible for the rf mass for the first 
time has been pointed out by 't Hooft [24]. Quantitatively this problem due to instantons 
was investigated in our previous work [25]. 

The formalism developed there will be generalized here in order to investigate this prob- 
lem (the large mass of the 7/ meson even in the chiral limit) quantititatively in quantum 
field theory as well in particular QCD within the above-mentioned zero momentum modes 
enhancement (ZMME) quantum model of the QCD nonperturbative vacuum [26,27]. 

However, a few general remarks in advance are in order. The quantum part of the VED 
is, in general, determined by the effective potential approach for composite operators [28] 
(see also Ref. [29]). It allows us to investigate the nonperturbative QCD vacuum, since in the 
absence of external sources the effective potential is nothing else but the VED. It gives the 
VED in the form of the loop expansion where the number of the vacuum loops (consisting of 
the confining quarks with dynamically generated quark masses and nonperturbative gluons 
properly regularized with the help of ghosts (if any)) is equal to the power of the Plank 
constant, Ti. As was underlined above, this quantity in general is badly divergent at least 
as fours power of the ultraviolet cutoff (for detail discussion see Shifman's contribution in 
Ref. [11] as well as our paper [30]). This reflects the fact (situation) that real QCD vacuum 
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being nonperturbative, nevertheless contains excitations and fluctuations of the gluon field 
confugurations there which are of pure perturbative character and magnitude. In order to 
deal with its true nonperturbative structure, all kind of perturbative contributions should 
be subtracted from the VED. In this way one obtains the truly nonperturbative VED which 
is precisely the one of main characteristics of the QCD true ground state in continuum 
theory. This is absolutely similar to lattice approach where by using different "smoothing" 
techniques such as "cooling" [31], "cycling" [32], etc. it is possible to "wash out" all types 
of the perturbative fluctuations and excitations of the gluon field configurations from the 
QCD vacuum in order to deal only with its true nonperturbative structure, i.e., free from 
all kinds of perturbative "contamination" ("noise"). 

We have already formulated [14,15] a general method how to correctly (in a manifestly 
gauge-invariant way) calculate the truly nonperturbative VED in the QCD quantum models 
of its ground state using the above-mentioned effective potential approach for composite 
operators [28]. The truly nonperturbative VED was defined as integrated out of the truly 
nonperturbative part of the full gluon propagator over the deep infrared (IR) region (soft 
momentum region). The nontrivial minimization procedure which can be done only by the 
two different ways (leading however to the same numerical value (if any) of the truly nonper- 
turbative VED) makes it possible to determine the soft cutoff in terms of the corresponding 
nonperturbative scale parameter which is inevitably present in any nonperturbative model 
for the full gluon propagator.^] The analysis of the truly nonperturbative Yang-Mills (YM) 
VED after the scale factorization provides an exact criterion for the separation of " stable vs. 
unstable" vacuum models. If the chosen Ansatz for the full gluon propagator is a realistic 
one, then our method uniquely determines the truly nonperturbative YM VED, which is 
always finite, automatically negative and it has no imaginary part (stable vacuum). In the 
same way (apart from some details concerning the chiral limit physics) can be determined 
the contribution to the truly nonperturbative VED from confining quark degrees of freedom 
[26,27]). 

The above-briefly-described general method can serve as a test of QCD vacuum different 
not only quantum, lattice [14,33] but classical models [15,34] as well. By applying our method 
to the classical theory, we thereby investigating the stability vs. instability of the vacuum of 
the classical models against quantum corrections. However, it is worth emphasizing that the 
only way to calculate the truly nonperturbative VED in quantum field theory in particular 
QCD from first principles [14,15] is the effective potential approach for composite operators 
[28] by substituting there a well-justified ansatz for the full gluon propagator since an exact 
solution(s) to its Schwinger-Dyson (SD) equation is not yet known. Moreover, it seems to 
us that it cannot be found in principle because of too complicated mathematical structure 
of the corresponding SD equation [2]. Fortunately, however, in order to calculate the truly 
nonperturbative VED it suffices to know its deep IR asymptotics which for realistic models 
usually coincides with its truly nonperturbative part [14]. 

In our previous works [14,26,27] we have formulated a new, quantum model of the QCD 
ground state (its nonperturbative vacuum): the ZMME model or simply ZME since we al- 



1 Any deviation of the full gluon propagator from the free perturbative one automatically assumes 
its dependence on a characteristic mass scale parameter [30]. 



4 



ways work in the momentum space. It is based on the existence and importance of such kind 
of the nonperturbative, topologically nontrivial quantum excitations of the gluon field con- 
figurations (due to self- interaction of massless gluons only, i.e., without involving any extra 
degrees of freedom) which can be effectively, correctly described by the g _4 -type behaviour 
of the full gluon propagator in the deep IR domain. Such strong IR singular behavior of the 
full gluon propagator can be reffered to as the strong coupling regime, the so-called "infrared 
slavery" [2]. In realistic models for the full gluon propagator its truly nonperturbative part 
usually coincides with its deep IR asyptotics, emphasizing thus the strong intrinsic influence 
of the IR properties of the theory on its nonperturbative dynamics. By applying the above- 
described general method to this model, the nonperturbative chiral QCD vacuum was found 
stable, i.e., having a "stationary" state, which can be a manifestation of a possible existence 
of a ground state in this model. Consequently, the corresponding truly nonperturbative YM 
VED is finite and negative. 

Within the ZME quantum model of the QCD ground state, the truly nonperturbative 
VED depends on a scale at which the nonperturbative effects become important. If QCD 
itself is confining theory, such a characteristic scale should certainly exist. The confining 
quark part of the VED depends in addition on the constant of integration of the correspond- 
ing quark SD equation in the chiral limit. The numerical value of the nonperturbative scale 
as well as the above-mentioned constant of integration is obtained from the bounds (1.4) 
for the pion decay constant in the chiral limit by implementing a physically well-motivated 
scale-setting scheme [26,27] .0 We have obtained the following numerical results for the truly 
nonperturbative VED due to ZME quantum model: 



where, obviously, the first and second values are due to upper and lower bounds in Eq. (1.4), 
respectively. Here and further on below Nf is the number of light flavors. We see that the 
confining quark part is approximately one order of magnitude less than the YM part and it 
is of opposite sign. 

It is instructive to compare our values (2.1) with those which are due to instantons in 
the chiral limit [25]. Within the above-mentioned RILM, for dilute ensemble, the truly 
nonperturbative VED is ej = —(6/4) x no- Here b = 11 — (2/3)iV/ is the first coefficient 
of the (3 function (see below) and uq is the instanton number density in the chiral limit 
[23,25]. Due to all reasonable estimates of this quantity (which follows from phenomenology 
or lattice approach) its numerical value cannot exceed its phenomenological value, i.e., Uq < 
n = 1.0 /m -4 . Thus, at maximum the truly nonperturbative VED due to RILM is 



2 In nonperturbative QCD all numerical results depend on the numerical value of the nonpertur- 
bative (characteristic) scale and not on a s as in phenomenological and perturbative QCD. This is 
due to phenomenon of "dimensional transmutation" in field theories with spontaneous symmetry 
breaking [35,36] (see also discussion in Ref. [30]). 



e ZME = -(0.01413 - JV/0.00196) GeV A } 
e ZME = -(0.009 - V/0.00124) Gel/ 4 , 



(2.1) 



ej = -(0.004179 - JV/0.00025) GeV 4 . 



(2.2) 
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It is clear that our values (2.1) approximately one order of magnitude bigger than those of 
instantons can provide at all (2.2). 

One of the main purposes in this paper is to generalize a formalism, developed earlier in 
our paper [25] , in order to directly calculate the gluon condensate as a function of the truly 
nonperturbative VED, the topological susceptibility and the mass of the r( meson in the 
chiral limit as a functions of the truly nonperturbative YM VED. It is based on using the 
trace anomaly (section 3) and low energy (section 4) relations, as well as the WV formula 
for the mass of the rj' meson in the chiral limit (section 5). It was directly applied to the 
ZMME quantum model of the QCD nonpertubative vacuum. In section 6 we present our 
discussion and conclusions. Our numerical results are shown in Tables I-III. 



III. THE TRACE ANOMALY RELATION 

The truly nonperturbative VED is important in its own right as one of the main char- 
acteristics of the QCD nonperturbative vacuum. Furthermore it assists in estimating such 
an important phenomenological parameter as the gluon condensate, introduced in the QCD 
sum rules approach to resonance physics [37]. The famous trace anomaly relation [16] in the 
general case (nonzero current quark masses m") is 

© W = ^G%G; U + £ m%q f . (3.1) 

where 9 W is the trace of the energy-momentum tensor and G a being the gluon field strength 
tensor while a s = g 2 /Ait. Sandwiching Eq. (3.1) between vacuum states and on account of 
the obvious relation (0|O MAt |0) = 4e t , one obtains 

4e t = (0|^G^GyO) +J2m o f (0\q f q f \0), (3.2) 

where (0|gjg/|0) is the quark condensate and e t is the sum of all possible independent, truly 
nonperturbative contributions to the VED (the total VED). In quantum theory it consists of 
two parts: the truly nonperturbative gluon part properly regularized with the help of ghosts 
(if any), e g , the so-called YM part at Nf = 0, and confining quark part with dynamically 
generated quark masses, e q , i.e., 

e t = t g + N f e q . (3.3) 

In general case it is impossible to use the weak coupling limit solution to the f3 finction, 
so it is convenient to introduce the gluon condensate as follows: 

G 2 ^ (G 2 ) = -(^G 2 ) = -(OI^G^JO). (3.4) 

This is nothing else but the most general parametrization of the gluon correlation function, 
(OIG^G^JO) (for details, see sect. VI below) which comes from the trace anomaly ralation 
(3.2) itself. If confinement happens then the /3 function is always in the domain of attra- 
tion (i.e., always negative) without IR stable fixed point [2]. Thus the gluon condensate, 
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parametrized as in Eq. (3.4), is always positive as it should be. The trace anomaly relation 
(3.2) then can be written as follows: 

4ei = -G 2 + $>°(0|g /g/ |0>. (3.5) 

/ 

Since the gluon condensate G2 is always positive and finite while quark condensate is always 
negative and also finite, then the truly nonperturbative total VED e t is also always negative 
and finite without having imaginary part since G2 and quark condensate are real finite 
numbers as it has been emphasized above. As was underlined in the preceding section, the 
confining quark part of the total vacumm energy density e q (which is always positive) as usual 
an order of magnitude less than the corresponding YM part e g (which is always negative), 
so the total truly nonperturbative VED e t is always negative in complete agreement with 
Eq. (3.5). On account of Eq. (3.3), it gives the gluon condensate (3.4) as a function of Nf 
if somebody knows, of course, how to calculate et and quark condensate from first principles 
(see below and our papers [14,26,27]). 

In the chiral limit (mj = 0) things drastically simplify. From Eq. (3.5) one obtains, 
= — 4e°, indeed, where superscript "0" means the chiral limit. Thus in this limit the 
truly nonperturbative total VED is nothing else but the gluon condensate apart from the 
overall numerical factor. In what follows we will saturate the total VED in this equation by 
our values (2.1), i.e., to put e° = €zme + ■■■■ The numerical results are shown in Table I. 



IV. THE TOPOLOGICAL SUSCEPTIBILITY 

One of the main characteristics of the QCD nonperturbative vacuum is the topological 
density operator (topological susceptibility) in gluodynamics (Nf = 0) [3] 

Xt = Kmi J d A xe^j^(0\T{q(x)q(0)}\0), (4.1) 

where q(x) is the topological charge density, defined as q(x) = (a s /47r)F(x)F(x) = 
(a s /An)F^ I/ (x)F^ u (x) and F^ v (x) = (l/2)e^ /pa Fp a (x) is the dual gluon field strength tensor 
while N c is the number of different colors. In the definition of the topological susceptibil- 
ity (4.1) it is assumed that the corresponding regularization and subtraction of all types 
of the perturbative contributions have been already done in order Eq. (4.1) to stand for 
the renormalized, finite and the truly nonperturbative topological susceptibility (see Refs. 
[3,17,18,38]). The anomaly equation in the WV notations is 

= N f (2/N c )(a s /4n)FF. (4.2) 

The topological susceptibility can be related to the nonperturbative gluon condensate via 
the low energy "theorem" in gluodynamics proposed by NSVZ [17] (by using the dominance 
of self-dual fields hypothesis in the YM vacuum) as follows: 

Kmi / *xe*(0\T{^GG(x)^GG(0)}\0) = -^(^G 2 ). (4.3) 
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Quite recently it was discussed by HZ in Ref. [18] (see also references therein) who noticed 
that it is not precisely a Ward identity, but rather is a relation between the corresponding 
correlation functions, indeed. That is why in what follows we call Eq. (4.3) as low energy 
relation or NSVZ-HZ relation. Thus there exist two proposals how to fix the numerical 
value of the coefficient £. The value £ = 2/6, where here and everywhere below (apart from 
section 6) 6 = 11, was suggested by NSVZ [17] who used the above-mentioned dominance 
of self-dual fields hypothesis in the YM vacuum. A second one , £ = 4/36, was advocated 
very recently by HZ using a one-loop connection between the conformal and axial anomalies 
in the theory with auxiliarly heavy fermions [18] (and references therein). For completeness 
we will use both values for the £ parameter]^ Let us note only that there exists an obvious 
relation between these two values, namely = (2/3)^nsvz- 
The anomaly equation in the NSVZ-HZ notations is 

<V5 = N f (a s /4n)GG, (4.4) 

with Nf = 3. Thus in order to get the topological susceptibility in the WV form from 
the relation (4.3), it is necessary to make a replacement in its left hand side as follows: 
GG — > (2/N c )FF in accordance with the anomaly equations (4.2) and (4.4). Then the 
WV topological susceptibility (4.1) finally becomes 

Xt = -£ 2 (^G 2 ) = -(2£)Vm, (4.5) 

where the second equality comes from Eqs. (3.4) and (3.6) by denoting the truly nonpertur- 
bative VED at Nf = as eym ■ The significance of this formula is that it gives the topological 
susceptibility as a function of the truly nonperturbative VED for pure gluodynamics, eyu- 
For numerical results in the NSVZ mode see Table II. In order to obtain numerical results for 
the topological susceptibility in the HZ mode, it suffices in this Table to multiply numbers 
in GeV units by factor (2/3) 2 and numbers in MeV units by factor y^2/3. 

V. THE 17(1) PROBLEM 

The topological susceptibility (4.1) assists in the resolution of the above-mentioned U(l) 
problem [1,2] via the WV formula for the mass of the rf meson (1.1). Within our notations it 
is expressed (in the chiral limit) as follows: f^rnl, = ANfXt, where f v i is the rj' residue defined 

in general as (0| E,=m,» QWqW) = *^N~ffv'P» and (Q\N f ^FF\rf) = (N c ^N~ f / 2)^171^ 
[3]. Using the normalization relation f v ' = y/2F®, one finally obtains 



3 In the weak coupling limit solution to the function (see, expression (6.2) below), the NSVZ-HZ 
low energy relation (4.3) for the NSVZ value of the £ parameter is reduced (as it should be) to that 
which was used in our previous paper [25] where the instanton vacuum in the chiral limit in various 
modifications was investigated. At the same time, the functional dependence of the topological 
susceptibility on the vacuum energy density is not, of course, changed (compare expression (3.5) 
of Ref. [25] with Eq. (4.5) below, on account of the corresponding values of the £ parameter). 
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Flml, = 2N fXt . 



(5.1) 



Eq. (4.5) then implies 

m%, = -2N f (§-) 2 e YM , (5.2) 

which expresses the mass of the rj' meson as a function of the truly nonperturbative YM 
VED. In previous expressions we omit for simplicity the superscript "0" in the pion decay 
constant as well as in m 2 ,. Our numerical results for the mass of the rf meson in the chiral 
limit (5.2) are shown in Table III. 



VI. DISCUSSION AND CONCLUSIONS 

A. Nf dependence of the gluon condensate 

Let us emphasize, that the general parametrization of the gluon condensate, introduced 
in Eq. (3.4), can be formally applied to any (3 function: to its weak coupling limit solution 
which we certainly know (see below) or to its strong coupling limit one which we certainly 
do not know yet. Also it is remains relevant whether one considers the chiral limit case or 
not. In phenomenology, however, the oftenly used parametrization of the gluon condensate 
is [37] 

G 2 ^ (G 2 ) = {—G 2 } = (OI^G^JO), (6.1) 

which in what follows we will call standard (or phenomelogical) parametrization of the gluon 
condensate. Of course, the gluon condensate in principle can be parametrized by different 
ways. Here it makes sense to remind that by parametrization it is understood which numer- 
ical factor (a s /ir, a s , etc.) is chosen to be associated with the gluon correlation function 
(0|G^G^|0). By comparing these two parametrizations of the same gluon correlation func- 
tion, it becomes clear that the general parametrization G 2 , given in Eq. (3.4), is useful in 
the nonperturbative calculations from first principles as it has been explained above. At the 
same time, it is useless in phenomenology since we do not know the (3 function there. How- 
ever, in phenomenology in many important cases it is legitimated to use the weak coupling 
limit solution to the (3 function, 

a 2 2 
(3(a s ) = -b^- + 0(a 3 s ), b = U--N f , (6.2) 

for example in instanton calculus [23,25]. Then the general parametrization (3.4) is reduced 
to the phenomenological one, Eq. (6.1), as follows: 

G 2 = \xG 2 = (1.375 - 0.083iV7) x G 2 . (6.3) 
8 

This relation makes it possible to resolve the old-standing problem in QCD phenomenology. 
It is well known that the phenomenologically extracted value of the gluon condensate does 
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not explicitly depend on Nf. It is usually nonexplicitly assumed that it is relevant for all Nf 
by maching from heavy to light quarks [37]. Though in general it should be Nf explicitly 
dependent quantity as it follows from the trace anomaly relation (3.5) on account of Eq. 
(3.3). This is true in the weak coupling regime and in the chiral limit as well. 

Our proposal is to consider (to interpret) G 2 as a real phenomenological condensate 
which on account of any numerical value of the standard gluon condensate G 2 , through the 
relation (6.3), becomes explicitly dependent on Nf, indeed. In other words, our proposal 
allows one to restore the explicit Nf dependence of the phenomenologically determined 
gluon condensate. The relation (6.3) in fact prescribes which numerical factor, explicitly 
depending on Nf, should be multiplied on G 2 in order to fix the correct Nf dependence of 
the gluon condensate in phenomenology. This interpretation is in agreement with the trace 
anomaly relation, of course. 

Another problem is the numerical value of the standard gluon condensate G2 itself. It can 
be taken either from phenomenology or lattice simulations (see below). The phenomenologi- 
cal analysis of QCD sum rules for the standard gluon condensate implies, G 2 ^ 0.012 GeV 4 , 
which can be changed within a factor of two [37]. However, it has been already pointed out 
[39] that QCD sum rules substantially underestimate the value of the standard gluon con- 
densate. The most recent phenomenological calculation of the standard gluon condensate 
is given by Narison in Ref. [40], where a brief review of many previous calculations is also 
presented. His analysis leads to the update average value as G 2 = (0.0226 ± 0.0029) GeV 4 . 
In Ref. [41] from the families of J/*f> and T mesons substantially larger values were recently 
derived, namely 0.04 < G 2 < 0.105 (GeV A ). Comparable with these bounds, the values for 
the standard gluon condensate were reported in recent lattice simulations [42]. 

Obviously that the standard gluon condensate in the chiral limit cannot be equal to its 
any phenomenological (empirical) value. Apparently it should be less, i.e., G® ~ v x G 2 , 
where v < 1 is some real number. Then the relation (6.3) in the chiral limit becomes 
G° 2 = (b/8) x G° 2 ~ (1/6/8) x G 2 = (1.375 - 0.083A^/) x v x G 2 . In Ref. [17] it has been 
argued indeed that the gluon condensate in the chiral limit is approximately two times 
less than its any phenomenological (empirical) value, G 2 ~ O.5G2, i.e., v ~ 0.5. Then the 
previous relation becomes G\ ~ (b/lQ) x G 2 = (0.6875 — 0.0415iVj) x G 2 , whatever numerical 
value of G 2 is. 

Our values for the gluon condensate (3.4) in the chiral limit (Table I) substantially larger 
than it is possible to estimate from phenomenology at all with the help of the above-derived 
relations. This difference is not only due to different physical observables as was noticed in 
Ref. [23]. It seems to us that it reflects rather different underlying physics as well. Our gluon 
condensate (which was calculated from first principles) is the strong coupling regime result 
and reflects the nontrivial topology of the true QCD vacuum where quantum excitations of 
the gluon fields play an important role. Precisely these types of gluon field configurations 
are mainly responsible for quark confinement and dynamical breakdown of chiral symmetry 
[14,26,27]. At the same time, phenomenologically extracted gluon condensate being the 
weak coupling limit result can be associated with classical instanton-type fluctuations in the 
QCD vacuum which by themselves do not confine quarks [32,43,44]. 

Concluding this subsection, let us note that the confining quark condensate contribution 
to the trace anomaly relation (3.5) vanishes in the chiral limit. However, due to all reasonable 
estimates of light quark masses, numerically its contribution is at 20% and thus comparable 
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with the systematic error in any phenomenological determination of the standard gluon 
condensate itself [37,45]. 

B. Some general remarks on the trace anomaly relation 

It is widely believed that the truly nonperturbative VED is determined by the trace 
anomaly relation (3.2) or equivalently (3.5). This may be so in phenomenology but not, of 
course, in nonperturbative QCD, where the situation is completely opposite (see discussion 
below). If by some phenomenological analysis it is possible to extract numerical values of 
the gluon and quark condensates then substituting them into the right hand side of the 
trace anomaly relation (3.5), one is able to estimate the total VED, indeed. However, 
this number such obtained tells a little (or even nothing) about detail structure of the QCD 
vacuum. It is not surprising since the gluon and quark condensates being its average (global) 
characteristics cannot account for its microscopic, detail structure. 

The quantity which is responsible for its detail structure is precisely the total truly 
nonperturbative VED which, by definition, is the sum of all possible, independent contribu- 
tions. In this way it refrects the complexity and veriety of the quantum, dynamical degrees 
of freedom (quarks and gluons) in the QCD vacuum. It may cotain contributions even from 
classical field configurations (instantons, for example). Precisely the detail structure of the 
QCD vacuum determines its global characteristics, and not vice versa. We have already 
explained in section 2 how the truly nonperturbative VED should be calculated from first 
principles in nonperturbative QCD. It is well known how to calculate from first principles 
the quark condensate as integrated out of the trace of the truly nonperturbative quark prop- 
agator over the deep IR region. Then the trace anomaly relation (3.5) becomes simply an 
algebraic equation for calculating the gluon condensate and in the chiral limit it is nothing 
else but the truly nonperturbative VED itself. 

The contribution to the right hand side of the trace anomaly relation (3.5) from the 
confining quark condensate vanishes in the chiral limit. However, it does not mean that the 
contribution from confining quarks with dynamically generated quark masses also vanishes 
in the left hand side of the trace anomaly relation. Even in the chiral limit the vacuum 
quark loops provide nonzero contribution into the total VED. Thus this quantity contains 
much more information about QCD vacuum than its global chracteristics can provide at 
all. Reflecting its detail structure the total VED is the main characteristics of the QCD 
nonperturbative vacuum and thereby is responsible for its global characteristics as well. 

C. Summary 

Using the trace anomaly relation (3.2), The NSVZ-HZ low energy relation (4.3) and 
WV formula for the mass of the v( meson (5.2) in the chiral limit, we have developed a well- 
justified formalism in the most general form in order to express the topological susceptibility, 
the gluon condensate and the mass of the r/ meson in the chiral limt as a functions of the 
truly nonperturbative VED. The crucial role in this belongs to the NSVZ-HZ low-energy 
relation (4.3). It allows one to relate the important topological quantities such as the gluon 
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condensate, the truly nonperturbative VED, the topological susceptibility, etc. to each other 
in a self-consistent way by providing well-justified coefficients between them. 

This formalism has been immediately applied to the investigation of the chiral QCD 
vacuum structure within its quantum ZMME model. Our numbers are collected in Tables 
I-III. Our values for the topological susceptibility in the NSVZ mode (Table II) slightly 
overestimate its phenomelogical value while in the HZ mode underestimate it since recalling 
£hz — {2/3)£,nsvz (see also text at the very end of section 4). The same situation takes 
place with the mass of the rf meson in the chiral limit (Table III). However, it is worth 
emphasizing that the precise validity neither of the WV formula (1.1) nor the NSVZ-HZ 
low energy relation (4.3) is not known. That is why we cannot address the question of 
systematic error bars one might assign to the final numerical values summarized in Tables II 
and III. The topological susceptibility to leading order in the large N c limit is of order N®. 
The next-to-leading correction of order TV" 1 , however maight be not very small at iV c = 3 
which was used in our calculations. 

Any way, the ZMME values of the truly nonperturbative VED (2.1) are of the necessary 
order of magnitude to account for the phenomenological value of the topological suscep- 
tibility and therefore to saturate the large mass of the rj' meson. In other words, the IR 
singularities likely to be presented in the QCD nonperturbative vacuum and summarized by 
the corresponding behavior of the full gluon propagator in the deep IR domain are mainly 
responsible for the large mass of the 7/ meson and consequently for the phenomenological 
value of the topological susceptibility. Precisely this type of gluon field configurations is 
closely related to quark confinement as well [14,26,27]. Also our proposal how to restore 
the correct Nf dependence of the phenomenologically extracted value of the standard gluon 
condensate seems to be important in general and for instanton calculus [23] in particular. 
It is well known that the gluon condensate due to instantons does not explicitly depend on 
Nf [23,25,37]. 

The author is grateful to H. Toki for many interesting and useful discussions on these 
topics during his stay at RCNP. It is also a pleasure to thank Gy. Kluge for many useful 
remarks and help. 
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TABLES 

TABLE I. ZMME model values for the gluon condensate. 



F° (MeV) 


G% (GeV 4 ) N f = N f = l N f = 2 


N f = 3 


93.3 
83.2 


0.05652 0.04936 0.04084 
0.0360 0.03104 0.02608 


0.0330 
0.02112 


TABLE II. ZMME model values for the topological susceptibility. 


F° (MeV) 


Xt (GeV 4 ) 


xT (MeV) 


93.3 
83.2 


0.00187 
0.00119 


207.9 
185.7 


TABLE III. 


ZMME model values for the mass of the rf meson in the chiral limit (MeV units) . 




m° (NSVZ) 


m°, (HZ) 


93.3 
83.2 


1135 
1015 


756.86 
677 
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